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The lattice thermal conductivity of crystalline Si nanowires is calculated. The calculation uses 
complete phonon dispersions, and does not require any externally imposed frequency cutoffs. No 
adjustment to nanowire thermal conductivity measurements is required. Good agreement with 
experimental results for nanowires wider than 35 nm is obtained. A formulation in terms of the 
transmission function is given. Also, the use of a simpler, nondispersive "Callaway formula", is 
discussed from the complete dispersions perspective. 
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Determining the thermal conductivity of semiconduct- 
ing nanowires plays a crucial role in the development of 
a new generation of thermoelectric materials^. However, 
experimental results are still scarce^. The first measure- 
ments of the thermal conductivity of Si nanowires in the 
~ 20 — 115nm range have been recently reported in ref|J 
From a theoretical point of view, it is important to be 
able to quantitatively calculate lattice thermal conduc- 
tivities of nanowires in a predictive fashion. In the past, 
several techniques have become widespread in the calcu- 
lation of thermal conductivities of bulk materialsa^iSii. 
These works are based on linearized dispersion models 
that involve a certain number of adjustable parameters. 
In principle, one might hope that a proper account of the 
boundary scattering in the nanowires might enable us to 
predict nanowire thermal conductivities using the proce- 
dures of ref. or ref. 0, without any parameters other 
than those obtained from fitting the bulk data. Never- 
theless, calculations using those two approaches did not 
yield reasonable results for nanowires^. A first goal of 
this paper is to show that, using accurate phonon dis- 
persion relations, predictive calculations of the nanowire 
thermal conductivities can be done, which do not involve 
any adjustable parameter to fit the nanowire experimen- 
tal curves. A second result in this paper shows that the 
use of a " modified" Callaway approach is still possible for 
nanowires, although this second approach is not predic- 
tive, and requires adjustment to nanowire measurements. 

For a nanowire suspended between two thermal reser- 
voirs, its thermal conductance, <j(T), is defined as the 
flow of heat current through the wire when a small tem- 
perature difference AT exists between the reservoirs, di- 
vided by AT and taking the limit AT -> 0. If the 
phonons transmit ballistically, the conductance is 
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where a is a set of discrete quantum numbers labeling 
the particular subbands in the one dimensional phonon 
dispersion relations, fs is the Bosc distribution, and v z 
is the speed of the phonon in the axial direction of the 



nanowire, z. Using v z dk z — dui we have^ 

where u4 is the lower (upper) frequency limit of sub- 
band a. Here we have defined the ballistic transmission 
function, S(,(a;), which simply corresponds to the number 
of phonon subbands crossing frequency ui. 

However, in general, phonon scattering takes place in 
the wire and at the contacts, so that the transmission 
function is not given by the simple mode count stated 
above. In general, the transmission function can be ob- 
tained from the Green function of the system. This al- 
lows to study the interesting issues of ballistic versus 
diffusive transport^, or the effect of one dimensional 
localization. 1 - We shall not be concerned with those phe- 
nomena here, but shall assume that diffusive transport 
takes place and the Boltzmann transport equation is 
valid. From this standpoint, it has been shown that a 
mode's transmission can be described in terms of a char- 
acteristic relaxation length in the wire's direction, \ a {u>), 
such that a diff (T) ~ £ a / h^A^A ^ Vz (a, k z )dk z 
where L is the nanowirc's length 12 . The thermal conduc- 
tivity is defined as fc(T) = —a, where s is the cross section 
of the nanowire. Therefore, in this diffusive regime 
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To obtain the thermal conductivity we must compute 
the complete dispersion relations for the wire, ui a (k z ), 
which, using v z (a, k z ) = d "2it* > , allows us to express 
the previous equation as 



with 



v ' 0, otherwise. 



If the modes' lifetimes r(w) are independent of the 
subband, and only depend on the frequency, then, the 
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relaxation lengths are given by 

X a (u)) = v z (ol,u)t(u)), uj l a <uj<uj f a (6) 

Knowing u> a (k z ) and t(oj), Eqs. i) ll(il) already allow us 
to compute the thermal conductivity. Nevertheless, it is 
useful to recast Eq. (@J into an equivalent form, so as to 
explicitly retain the " transmission function" concept. For 
this, we define N b (ui) as the number of phonon subbands 
crossing frequency 

N b (u) = e(w-wi)e(^-w), (?) 

a 

The average group velocity in the axial direction is de- 
fined as 

(v,{w)) = fcv t {a,w)J /N b (u) (9) 

where the sum extends to subbands a such that cj* < 
lo < uj^. From Eqs. (|5I6|) and 10, it follows that 
E^«H = T ( UJ )J2 a v A a , UJ ) = r(uj)N b (uj)(v z (uj)}. 
Thus Eq. can be recast into its completely equiva- 
lent formSi 

k{t) -sL SM ^ dw (l0) 

with — H(w) = N b (ui)T(ui)(vJuj)). (11) 
s 

Here, the diffusive transmission function 5 has been de- 
fined, and 

N b (u) = N b (cj)/s. (12) 

Computation of the transmission function requires the 
obtention of the full phonon dispersions for the nanowire. 
We now proceed to separately study each of the three 
factors in Eq. (|11|) : N b (uj), t(lu), and (v z (w)). 

N b (uj) is obtained from the complete dispersion rela- 
tions for the wire. To compute them, we have used the 
interatomic potential proposed by ref. Il3l in which the 
energy of the system is given as the sum of two and three 
body terms. The two body terms are defined as 

5E Q (i,j) = | C (dij ~d ) 2 /d 2 

for every pair of nearest neighbors i and j, where dij 
is the distance between the atoms and do is the lattice 
equilibrium distance. The three body terms are defined 
as 

6E 1 (i,j,k) = \C 1 8&l j , k 

for every pair of bonds joining atoms i, j and k, where 
<50 is the deviation with respect to the equilibrium angle 
between the two bonds in the lattice. The constants are 



ion 




0.2 0.4 0.6 0.8 1 

co [10 2 THz] 



FIG. 1: Solid line: function iV b °°(w) for the (110) direc- 
tion, calculated using complete phonon dispersions. Dotted 
line: parabolic approximation to N§°(uj), for a non dispersive 
medium. Inset: dispersion relations for a 2.2 nm wide wire 
with frozen boundary. 



C = 49. W and C x = 1.07eVii£. With this potential 
the dynamical matrix of the system is constructed, and 
the dispersion relations for the wire are calculated. The 
inset of fig. H shows one example. The bulk dispersions 
obtained with the Harrison potential are close to those 
experimentally measured. Although a more sophisticated 
calculation of the interatomic potential is possible using 
ab initio techniques, this would only result in minor dif- 
ferences in the calculated N b and the conductivities com- 
puted from it. 

The dispersions in general depend on the cross sec- 
tional shape of the wire, and on whether its surface is 
reconstructed or not, clean or coated by an overlayer, 
etc. As the surface/ volume ratio decreases, these differ- 
ences become smaller, and for wide enough wires of cross 
section s, the function N b (uj) = N b (u>)/s approaches a 
limiting form, A^°(o;), independent of the surface fea- 
tures. To ascertain whether it is valid to approximate 
N b (oj) by A^^°(ti>) for nanowires wider than 35 nm, the 
former was calculated for nanowires of increasing width, 
with frozen boundary conditions. (As compared to "free 
boundary" , frozen boundary conditions are more simi- 
lar to the experimental situation, where 1 nm of S1O2 
coats the wire.) For a romboidal cross section wire of 
17 nm side, N b deviates from the limiting value by 
less than 0.1%i2i Therefore, for the widths considered 
here (I > 35nm), N b = is a valid assumption. In the 
thermal conductivity calculations presented here, Nj^ for 
the 110 direction was used (fig-CJ. 

The phonon lifetime is commonly given by Matheis- 
sen's rule, expressing the total inverse lifetime as the sum 
of the inverse lifetimes corresponding to each scattering 
mechanism. For bulk or a macroscopically thick whisker, 
the expressions used for boundary, anharmonic and im- 
purity scattering arei^: 

r b - 1 = ({l/c) b lF)-i 
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T- l = BTu 2 e- c ' T (13) 

From here on, ()& denotes averaging in the three acous- 
tic branches, i.e. (l/c)b = 5 X)i=i l/ c ii where c, is the 
speed of sound for each branch. (Do not confuse with 
averages over subbands, denoted by ().) A, B and C are 
numerical constants specified later. I is the wire's lateral 
dimension, and F relates the boundary scattering rate 
to the shape and specularity of the sample's boundary. 
Eqs. (|13fl are always valid for wire thicknesses larger than 
a certain threshold, L^rj- On the other hand, if the width 
of the whisker is very small, such that L << L^r> , the fre- 
quency difference between consecutive phonon subbands 
gets large. This confinement modifies the inter-subband 
scattering, and the validity of Eqs. (|13fl breaks down^ A 
criterion to estimate L^d is: confinement is important if 
the energy spacings between consecutive subbands, HAuj, 
are larger or comparable to the thermal energy fc^T. 
The largest intersubband spacings are of the order of 
TiAuj ~ T"^ c ' Thus, for Eqs. ^| to be valid at tem- 
peratures 40°K or higher, we need I > lOnm. This is 
well fulfilled for nanowires in the 37-115 nm range, as 
the ones experimentally measured^. Therefore, for these 
widths, we use Eas. 1131 without modification. The good 
agreement with experimental results corroborates their 
validity. For sub 20 nm widths, confinement effects may 
be expected^ 

Anharmonic scattering should include the effects of 
umklapp as well as normal processes. In some works 
normal processes have been omitted in favor of umklapp 
scattering, without extensive discussion^^. In general, if 
resistive processes dominate, the normal relaxation time 
can be disregarded 1 ^. For the higher frequencies, umk- 
lapp scattering is proportional to w 2 , and it thus domi- 
nates over normal processes, which depend linearly on 
frequency^. At the lower frequencies, in the case of 
nanowires, boundary scattering dominates. This justi- 
fies the form of r a in Eq. ^1 which is the functional de- 
pendence for umklapp scattering at the higher frequen- 
cies. Constants B and C were adjusted to reproduce the 
bulk material's experimental thermal conductivity curve, 
yielding B = 1.73 x lCT 19 s/A' and C = 137.3AT^i Pa- 
rameter A — 1.32 x 10~ 45 s 3 is analytically determined 
from the isotope concentration, and it should not be 
adjusted 1 ^, so we maintain the value given in 

The third term in Eq. lfTT|l . (v z (lj)), is also calculated 
from the dispersion relations. In a frozen boundary wire, 
it is a matter of simple geometry to prove that at the 
lower part of the spectrum^ 

(v z ) ~ ^l/c) b /<l/ C 2 ), (14) 

where (l/c 2 )b = § X)i=i V c f • The relaxation times, 
Eq. I|13l) are quite rough, since they only depend on the 
lower spectrum speeds of sound. Thus, no additional ap- 
proximation is incurred if Eq. (|14|) is used for the third 
term in Eq. (|11|) . 
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FIG. 2: Thermal conductivities versus temperature calculated 
using the complete dispersions transmission function, for Fl = 
1.05 x 37nm (solid), 1.3 x 56nm (dotted) and 1.15 x 115nm 
(dashed). Dots: experimental results from ref. |3J. Inset: bulk. 

Fig. El shows the results of the calculation using the 
complete dispersions transmission function into Eq. HUH . 
Calculations for nanowires of diameters Fl — 38.85, 72.8 
and 132.25 nm are shown, together with experimental re- 
sults from ref.OJ Although the TEM measured diameters 
of the experimental nanowires were 37, 56 and 115 nm, 
variations in cross sectional shape, as well as thickness 
of the nanowire's oxide surrounding layer are possible, 
so that the cross sectional area is not accurately known. 
Assuming / equal to the TEM diameters, values can be 
assigned to F, which in the current case are F ~ 1.05, 
1.3 and 1.15. These values suggest that the boundary 
scattering is very diffusive, rather than specular. 

Let us see what would happen if, instead of the atom- 
istically calculated complete dispersions, a nondispersive 
medium were used. In such a case, for frozen boundary, 
function Ng°(u>) is analytical, having a parabolic depen- 
dence: 

(non dispersive) N b °° (w) ~ £?=i J^ 1 2nk ± -^ 

^ 2 (l/c 2 ) fc (15) 

Substituting Eq. JTHJl into Eqs. JTJJ and d), tne Call- 
away formula*^ is obtained. Since now the subbands 
extend to infinite frequency, an upper cutoff has to be im- 
posed, which was not necessary when using the complete 
dispersions, because for these ones the upper bands' lim- 
its are well defined. Traditionally, the Debye frequency, 
ujd (86 THz for Si), has been used as cutoff 4 . If we do 
this, the calculated thermal conductivities for nanowires 
are in large disagreement with experimental measure- 
ments. As fig. Eta) shows, no calculated curve agrees 
in shape with any experimental curve: the theoretical in- 
flection points are always too high. Comparison of the 
parabolic approximation to N§° (Eq. I15f) with the cor- 
rect N£°, in fig. 2] shows that the former grossly overes- 
timates the transmission of the high frequency phonons, 
if its cutoff is set at uir>- Nevertheless, if a lower cut- 
off luc is chosen, the " Callaway" and " correct" Njf (us) 
curves are more similar. Using toe, the shapes of the 
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FIG. 3: Thermal conductivity normalized by its value at 
320° A", for nanowires of different widths, calculated by the 
Callaway formula with cutoff at lod (86 THz) (a), and at cue 
(42 THz) (b). (Experimental results from ref. |3-) Inset: bulk 
thermal conductivity for the two cases (overlapping almost 
completely) and experimental curve from ref. |5j (dots). 

theoretical thermal conductivity curves agree well with 
the experimental ones (fig- Oft))). The values of the life- 
time parameters B and C are those that yield the correct 
bulk limit: 3.9 x 10~ 19 s/K and 223. 1°K for fig-Efa), and 
1.7 x 10~ 19 s/K and 151. 8°K for fig.EJb). 

For bulk, it is possible to obtain a good fit of experi- 
mental data using either ujq or u> d (see inset of fig. |2l , 
because the dominant umklapp scattering makes the role 



of high frequencies less important. For nanowires, bound- 
ary scattering dominates, and high frequencies play a 
more important role. Therefore, using the experimen- 
tal Debye frequency as a cutoff for the Callaway formula 
does not yield correct results for nanowires, although 
good results can be obtained with a lower, adjusted cut- 
off. This one must be determined by comparison with 
nanowire experimental results {non predictive approach). 
On the other hand, if the atomistically calculated, com- 
plete dispersion relations are used, good agreement with 
experimental results is obtained (fig-OJ), with no need for 
prior knowledge of any experimental measurements for 
nanowires. 

To conclude, it was shown that, by using complete, 
atomistically computed phonon dispersions, it is possi- 
ble to predictively calculate lattice thermal conductivity 
curves for nanowires, in good agreement with experi- 
ments. On the other hand, it was not possible to obtain 
correct results with the approximated Callaway formula. 
Still, good results could be obtained with a modified Call- 
away formula, although non predictively. The results of 
this paper are only expected to apply to " nanowhiskers" 
for which phonon confinement effects are unimportant. 
Si nanowires wider than ~ 35nm are within this cate- 
gory. 

I am indebted to D. Li, A. Majumdar and Ph. Kim for 
communication of their experimental data, and discus- 
sions, and to A. Balandin and Liu Yang for discussions. 
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which agrees with our ob- 



tained value (7 = Gruneisen parameter, /1 = shear modu- 
lus, Vo = volume per atom, and ujd = Debye frequency.) 
Constant C has the order of magnitude of the transverse 
phonons' frequency near the zone edge—*—. 
25 For this, use u> a (k z ) ~ a + k\ , a = {i, k±}, and 
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